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ABSTRACT 

For cylindrical differentially rotating plasmas, we study large-scale magnetic field gen¬ 
eration from finite amplitude non-axisymmetric perturbations by comparing numerical 
simulations with quasi-linear analytic theory. When initiated with a vertical magnetic 
field of either zero or finite net flux, our global cylindrical simulations exhibit the mag- 
netorotational instability (MRI) and large scale dynamo growth of radially alternating 
mean fields, averaged over height and azimuth. This dynamo growth is explained by 
our analytic calculations of a non-axisymmetric fluctuation-induced EMF that is sus¬ 
tained by azimuthal shear of the fluctuating fields. The standard effect” (shear of 
the mean field by differential rotation) is unimportant. For the MRI case, we express 
the large-scale dynamo field as a function of differential rotation. The resulting radially 
alternating large-scale fields may have implications for angular momentum transport 
in disks and corona. To connect with previous work on large scale dynamos with local 
linear shear and identify the minimum conditions needed for large scale field growth, 
we also solve our equations in local Cartesian coordinates. We find that large scale 
dynamo growth in a linear shear flow without rotation can be sustained by shear plus 
non-axisymmetric ffuctuations-even if not helical, a seemingly previously unidentified 
distinction. The linear shear flow dynamo emerges as a more restricted version of our 
more general new global cylindrical calculations. 


1 INTRODUCTION 

Astrophysical rotators, such as stars, galaxies, and accretion 
disks, commonly show evidence for contemporaneous pres¬ 
ence of disordered turbulence and magnetic fields ordered on 
spatial or temporal scales larger than those of the fluctua¬ 
tions. Explaining this circumstance has been a long standing 
challenge. In situ amplification of large-scale magnetic fields 
via some type of large-scale dynamo is likely but how these 
dynamos operate and saturate in each context remains an 
active subject of research (for reviews see (Brandenburg & 
Subramanian 2005; Blackman 2015)). How such fields grow 
given the presence of fluctuations, what are the best analy¬ 
sis methods, and what minimum ingredients for growth are 
needed (Vishniac & Brandenburg 1997; Brandenburg et al. 
2008; Yousef & et al. 2008; Heinemann et al. 2011; Herault 
et al. 2011; Squire & Bhattacharjee 2015) are topics of active 
investigation. 

Beyond stellar and galactic contexts, evidence for large- 
scale field growth is seen in magnetically dominated labora¬ 
tory plasmas (Ji et al. 1995; Cothran et al. 2009), and in local 
and global simulations (Brandenburg et al. 1995; Ebrahimi 
et al. 2009; Lesur & Ogilvie 2010; Davis et al. 2010; Simon 
et al. 2011; Guan & Gammie 2011; Sorathia et al. 2012; 
Suzuki & Inutsuka 2014) of the magnetorotational instabil¬ 
ity (MRI) (Velikhov 1959; Balbus & Hawley 1991). Large- 


scale fields in MRI flows have been associated with the suste¬ 
nance of MRI turbulence(Lesur & Ogilvie 2008; Davis et al. 
2010; Simon et al. 2011) are correlated with the convergence 
of Maxwell stress (Guan & Gammie 2011; Nauman & Black¬ 
man 2014), and can influence corona formation (Blackman 
& Pessah 2009). For a single MRI mode, large-scale mag¬ 
netic fields generated via an EMF can cause MRI satura¬ 
tion (Ebrahimi et al. 2009). In short, the large-scale dynamos 
of MRI-unstable systems are of interest both as phenomena 
on their own, and because they may be closely connected 
to angular momentum transport in accretion disks by local 
and nonlocal Maxwell stresses (Blackman & Nauman 2015). 
In addition to numerical simulations, flow-dominated lab¬ 
oratory experiments are also investigating the MRI MHD 
unstable systems in Taylor-Gouette flow geometry (Good¬ 
man & Ji 2002; Rudiger et al. 2003; Kageyama et al. 2004; 
Noguchi et al. 2002; Sisan et al. 2004; Stefani et al. 2007). 

The electromotive force (EMF) from correlated veloc¬ 
ity and magnetic field fluctuations is important to all large 
scale dynamo theories (Moffatt 1978). In general, correlated 
fluctuations in the EMF facilitate large-scale field amplifica¬ 
tion and the form that this takes for cylindrical MHD shear 
flows is the focus of the present paper. 

As the role of large-scale dynamos and large scale fields 
for MRI turbulence and angular momentum has become in- 
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creasingly recognized, there is a need for truly global strat¬ 
ified simulations in the long run to best compare to real 
astrophysics disks. But choices must always be made both 
due to limited computational resources and for isolating key 
physical processes that contribute to the global dynamics. 
The shearing box model has been the workhorse for sim¬ 
ulating MRI turbulence for this purpose- but it is a local 
model and has limitations associated with boundary condi¬ 
tions, box size and so is a limited model for real astrophysical 
disks (Regev & Umurhan 2008; Bodo et al. 2008; Blackman 
& Nauman 2015). The cylindrical model used here also has 
some complementary limitations due to its boundary con¬ 
ditions, but on the other hand provides solutions in a real 
global domain. Moreover, as unstratified shearing simula¬ 
tions go, the shearing box in itself might be thought of as a 
more restrictive approximation of a global unstratified cylin¬ 
der of the sort used here. 

In contrast to previous studies of local Cartesian shear¬ 
ing box addressing large scale field growth from nonhelical 
turbulence and linear shear, (Brandenburg 2005; Yousef & 
et al. 2008; Heinemann et al. 2011; Squire & Bhattachar- 
jee 2015), we focus on the more general large-scale dynamo 
from the combination of non-axisymmetric perturbations 
and global differentially rotating flows in a cylinder. We 
show how the combination of imposed non-axisymmetric 
fluctuations and differential rotation, or linear shear of the 
fluctuating field, is sufficient to source the electromotive 
force and generate a large scale magnetic field in this cylin¬ 
drical geometry. We present the complete quasilinear form 
of the EMF and show that it models favorably the results 
direct numerical simulations (DNS) of the MRI when the 
magnitude and growth rate of these initial fluctuations in 
the simulations are used as inputs to the quasi-linear sin¬ 
gle model analysis. The single mode analysis proves useful 
in showing explicitly that mode-mode coupling is not essen¬ 
tial for growth, and for identifying which terms in the EMF 
dominate. We also show that these conditions for large scale 
field growth do not depend on whether the shear profile is 
favorable or unfavorable to the MRI as long as there is a 
physically motivated source of fluctuations. 

To identify the minimum requirements for large scale 
growth, to connect with previous work, and to compare with 
the global cylindrical case, we carry out analogous calcula¬ 
tions in local Cartesian coordinates. Comparing cylindrical 
and local Cartesian models, we find that in each case the 
fluctuation-induced EMF has separate contributions that 
depend respectively on 1 ) non-uniformity of the radially 
sheared non-axisymmetric perturbations 2 ) the background 
differential rotation (cylindrical) or linear shear (Cartesian) 
3) the mean angular velocity. These three vertical EMF 
terms can separately generate a large-scale magnetic field. 
We discuss them in the context of our numerical simulations 
and the minimum requirements for large scale field growth. 
We find that non-axisymmetric fluctuations plus linear shear 
OR uniform rotation provide the two most minimal combi¬ 
nations needed for large scale field growth in local Cartesian 
limit. 

In Sec. 2, we present evidence of large-scale toroidal 
fields from global nonlinear MHD DNS of the MRI in a cylin¬ 
drical setup. We derive the general form of the EMF in the 
quasilinear approximation in cylindrical geometry in Sec. 3. 
There we also discuss the role of the EMF in field growth 


(Sec. 3 .2) and consider an example limiting case where ra¬ 
dial gradients of the fluctuations are ignored (Sec. 3.3). We 
compare the EMF expression from cylindrical analytics with 
numerical calculations in Sec. 4. We find that the terms with 
direct dependence on mean differential rotation contribute 
most to the dynamo seen in the simulations. This is ex¬ 
amined for a specific nonaxisymmetric MRI mode and it is 
shown that a toroidal large-scale field is directly generated 
through the vertical EMF as the result of the coupling of 
a small-scale fluctuations with the differential rotation. The 
traditional “fl effect” (i.e. growth of mean toroidal field from 
mean poloidal field by shear) (Moffatt 1978; Parker 1979) is 
unimportant when the initial mean field is purely vertical. 
Visualizations of the field lines from our nonlinear MRI sim¬ 
ulations of the cylinder are presented in Sec. 5 to highlight 
why non-axisymmetric perturbations are needed for large 
scale field growth. In Sec. 6 we repeat the quasilinear analy¬ 
sis in local Cartesian coordinates and derive general forms of 
the EMF in this geometry. We discuss the associated impli¬ 
cations for field growth. By analogy to the specific simplify¬ 
ing example in Sec. 3.3, we discuss the limiting case in which 
the X gradients of the fluctuations vanish (Sec. 6.2). Finally, 
in Sec. 7, we emphasize that large-scale fields can also be 
generated even if the rotation profiles would imply stabil¬ 
ity to the MRI, as long as there is some external supply of 
non-axisymmetric fluctuations. We conclude in Sec. 8 , and 
also present a Table summarizing the ingredients needed for 
dynamo action. 


2 DIRECT NUMERICAL SIMULATIONS IN A 
CYLINDER 

We begin with our results from global DNS MHD simula¬ 
tions of the MRI in cylindrical (r, (j), z) geometry using the 
DEBS (Schnack et al. 1987; Ebrahimi et al. 2009) initial- 
value code to solve the nonlinear, viscous and resistive MHD 

= -E = S'V X B - 77 J (1) 

= -SpV.VY + S3 ><-& +- S^VF{2) 

= -S'V • (PV) - S'(r - 1)PV • V (3) 

= -SV ■ (pV) (4) 

= V X A (5) 

= V X B (6) 

where the variables, p, P,V, B, J,, and F are the density, 
pressure, velocity, magnetic field, current, and ratio of the 
specific heats, respectively. We use the same normaliza¬ 
tion (Schnack et al. 1987; Ebrahimi et al. 2009; Ebrahimi 
& Bhattacharjee 2014), where time, radius and velocity are 
normalized to the outer radius a, the resistive diffusion time 
Tn = /pop, and the Alfven velocity Va = Po/v^PoPo, re¬ 
spectively. The dimensionless parameters, S = tuVa/o. and 
Pm, are the Lundquist number and the magnetic Prandtl 
number (the ratio of viscosity to resistivity), respectively, 
the initial state satisfies the equilibrium force balance con¬ 
dition = pV^/r, where /?o = 2poPo/Po is normalized 

to the axis value, and the initial pressure and density profiles 


equation! 
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dt 
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dP 
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Figure 1. The generation of large-scale toroidal magnetic field 
(averaged over vertical and toroidal directions) in the zero-net 
flux simulations in the r-t plane, (a) when all the Fourier modes 
are included and (b) only one non-axisymmetric mode m=l is 
evolved (Rm=3100, Pm=l, /3o = 10®, Vq/Va = 31). 


are assumed to be radially uniform and unstratified. Pres¬ 
sure and density are evolved, however, they remain fairly 
uniform during the computations. A no-slip boundary con¬ 
dition is used for the poloidal flow and flow fluctuations. The 
inner and outer radial boundaries are perfectly conducting 
so that the tangential electric field, the normal component of 
the magnetic field, and the normal component of the velocity 
vanish. The tangential component of the velocity is the ro¬ 
tational velocity of the wall. The azimuthal (0) and axial {z) 
boundaries are periodic. We assume a radial pressure gradi¬ 
ent balances the centrifugal force in equilibrium, but radial 
gravity and a radial pressure force are interchangeable for 
our incompressible, unstratified circumstance. The pressure 
gradient, rather than gravity, is what balances the centrifu¬ 
gal force in cylindrical laboratory experiments designed to 
test the MRI (Goodman & Ji (2002)). 

All variables are decomposed as /(r, rj), z, t) = 
E(m.fc) = U{r,t)) + f(r, (jj z,t), where 

(/) is the mean (m = k = 0) component, and / is the 
fluctuating component. Mean quantities (indicated by 
brackets (()) or overbars) are azimuthally and axially 
averaged, but remain dependent on radius (r). Equations 
(1-6) are then integrated forward in time using the DEBS 
code. The DEBS code uses a finite difference method with a 
staggered grid for radial discretization and pseudospectral 
method for azimuthal and vertical coordinates. In this 
decomposition, each mode satisfies a separate equation 
of the form = Lrr,,kfm.k + ,k') ,k' ^ 

where Lm.k is a linear operator that depends on /o,o(r,t), 
and j, j./ is a nonlinear term that represents the 

coupling of the mode (m,k) to all other modes (m ,k).{ 



t/"T r 


Figure 2. Total magnetic energies, = 1/2 J B'^dr^ and Wr = 
1/2/ B^dr^, and large-scale toroidal magnetic energy, {W^) = 
l/2/(S^)2dr3 , vs. time for net-zero flux 3-D MRI computations. 


This latter term is evaluated pseudospectrally.) The time 
advance is a combination of the leapfrog and semi-implicit 
methods (Schnack et al. 1987). 

We initiate simulations with a Keplerian flow (V/(r)) = 
Vor~^^^ and uniform magnetic held B = Bgz (with non-zero 
initial net-flux) or B = Bosin(27r(r —ri)/(r 2 —ri))/r 2 , (with 
zero-net-flux) where ri, r 2 are the inner and outer radii. 
Fully nonlinear simulations with all Fourier modes included 
(with radial, azimuthal and axial resolutions of nr=220, 
0 < m < 43 and —43 < n < 43) show that large-scale 
magnetic helds are generated [Figure 1(a)]. In all of our 
simulations, the initially weak vertical magnetic field, gets 
redistributed, amplified at inner radii and reduced at outer 
radii. Initially =0, but a toroidal large-scale (averaged in 
<j) and z) held grows via the correlation of non-axisymmetric 
MRI-induced fluctuations. The sustenance of total toroidal 
and radial magnetic held energies, as well as the large-scale 
toroidal magnetic energy during the computation for 
net-zero flux are shown in Fig. 2. The radially dependent 
large-scale field [Figure 1(a)] is also sustained in time. 

The code can also be used to compute the nonlin¬ 
ear evolution of a single mode evolution for which the 
initial conditions consist of an equilibrium {f{r)) plus a 
single mode /m,fc(r, perturbation. The initial 

amplitude /m,fc(r, 0) is a polynomial in r that satisfies 
the boundary conditions at r = ri and r = r 2 . The 
initial amplitudes of all other modes are set to zero. Only 
the mode (m, k) is then evolved; however, the m = 0, 
k = 0 component (the background) is allowed to evolve 
self-consistently. The evolution of the background profile 
/o,o can affect the evolution of the mode (m, k) and cause 
the mode to saturate. In a fully nonlinear computation, 
all modes are initialized with small random amplitude 
and are evolved in time, including the full nonlinear term 

(J^^,k.m'.k')- 

To facilitate the analytic investigation of the large-scale 
held generation, we carried out nonlinear single-mode 
(i.e. single value of m and k) simulation of m=l non- 
axisymmetric MRI (in which only one fluctuation mode and 
the mean fields self-consistently grow) for zero-net-flux and 
net-flux configurations. As seen in Fig. 1(b) at t ~ 0.02, as 
the m=l MRI mode amplitude approaches saturation, a 
large-scale toroidal field is also generated. 
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3 QUASILINEAR ANALYTIC 

CALCULATIONS OF EMFS IN A CYLINDER 

3.1 Deriving the general form of the EMFs 


where 

£-ZVr' 


m'ykzBo 

{m? + 


rVrV*' 

7^ + 


( 13 ) 


To identify the origin of large-scale magnetic field growth in 
the DNS simulations, we employ quasilinear analytical cal¬ 
culations for the single-mode case. Given initial fluctuations 
we calculate the fluctuation-induced EMF £ = (V x B) from 
linearized eigenfunctions. The EMF is the source of large 
scale field growth. All averaged correlations are presented 
in terms of the radial Lagrangian displacement of a fluid 
plasma element, (Frieman & Rotenberg 1960), and the mean 
quantities. We assume perturbed quantities of the form 

+ for 

a cylinder of outer radius a and height L = 2-Elk. In the 
presence of an equilibrium mean flow, self-adjointness of the 
linear stability problem is lost (Frieman & Rotenberg 1960). 
Thus, for nonaxisymmetric modes (nonzero m), the eigen¬ 
values 7 c = 7 -I- iWr and the eigenvectors ^(r) can be com¬ 
plex, where 7 and Wr are the growth rate and the oscillation 
frequency of the mode, respectively. To isolate the role of 
shear flow on the dynamo effect in the quasilinear analytical 
calculations below, we impose an initial, uniform Bq. 

For a single MRI Fourier mode, the cylindrical coor¬ 
dinate components of the linearized momentum equation 
in terms of the Lagrangian displacement vector ^ (Chan- 
drasekar 1961) are 


[ 7 ^ +U3\ + 2rQ.{r)Q. (r)]^c - 2'^n{r)£^^ ^ ~~dr 

( 7 ^ -I- -I- 2^Q.{r)£,r = ( 8 ) 

[f+Lol)(,. = -ik.X, (9) 


where 7 = 7 -!- iu>r — imil{r), X = P + B ■ Bq, w\ = klB^/p 
and B(r) = V^(r)lr is the angular velocity. By including 
this imposed mean flow in the dehnition of the Lagrangian 
displacement vector, the velocity fluctuations in an Eulerian 
frame are given by V = d$,/dt -I- V x (^ x V), with compo¬ 
nents Vr = % = 7C<A - and 14 = 7 ^ 2 - 

For small resistivity, the magnetic held perturbations 
can be directly related to the displacement via B = ik^Bo^. 
Using this along with incompressibility and Eqs.( 7-9), we 
can eliminate X and the azimuthal and vertical displace¬ 
ments can be written in terms of as 




1 

1 -I- rrP/(r'^k'l) 


' - 2 B(r )7 ^ 


im 



( 10 ) 


and 




mk^r 

{r^kl + m?) 


( - 2 n(r) 7 \ 






i 

rkz 


(1 _ 

\ -I- m? 




( 11 ) 


where the primes indicate radial derivatives. The quasilinear 
EMF components fz = (V xB)z = ^Re{V*B^ — V^Br) and 
= (Y X B),/, = ^Re{V*Br — V*BG} can now be written 
in terms of the radial velocity fluctuations 14 = 7 Cr, 

£z = (V X B)z = £zvr' + £-zQ.' + £zn^ (12) 


£zn' — — 


rnr^k^Bo 


fl'(r)fl(r)r|U,| 


(72 + fl"(r ))2 


£zn = 


(m? -\- k^r^) 
2yklBon{r)En 2 , 772 


G 2 


-(7+11 -u;l)|14| 


myfczRo 

(m2 -I- fc|r2) 


1141 


72 -I- B^(r) 


and 

£<!> = (V X B),^ = £^vr' + £<i>n' + 
where 

-yBo w?^Bo ] (rl414* ) 


£ (j)Vr' 


£4>n' — — 


r r(m 2 -|- k^r^) 

72 -I- B^(r) 


7 B 0 

Q' {r)Q,(r)r 


r r{m? -\- klr"^) 

_ (72 + H"(r ))2 _ 


(14) 

(15) 

(16) 

(17) 

(18) 


_ 2-fklrBoQ.mQ{r) ^ 2 ,2 n2M77i2 

~ ^2 \^A 'J ** )|Vv| 


-yBo m^-yBo 

IK-P 

r r{m? -\- k^r^) 

72 -I- D^(r) 


where 

= (m^-|-fc2r^)[7^-|-B^(r)][47^B^(r)-|-(7^-|-B^(r)-|-a;i)^] 

( 20 ) 

and £l{r) = mQ{r)—uJr. In an ideal MHD cylindrical plasma, 
Eq. (12) provides the complete quasilinear form of the ver¬ 
tical huctuation-induced EMF in terms of radial perturba¬ 
tions. The hrst term on the RHS, £zvr', depends on the 
non-uniformity of the radial displacement of the mode. The 
second term £zni, which depends on the differential rotation 
f 2 '(’') is sufficient to directly produce a nonzero huctuation- 
induced dynamo term. The free energy source appears 
in this term. The third term, £zn, shows the dependence of 
the vertical EMF on angular velocity. 

The linearized cylindrical solutions (7 and ^r) for non¬ 
axisymmetric flow-driven and MRI modes have been pre¬ 
viously examined (Bondeson et al. 1987; Ogilvie & Pringle 
1996; Keppens et al. 2002). Here we do not solve the eigen¬ 
value problem to find the G for nonaxisymmetric modes 
but (in Sec. 4,) extract the linearized solutions directly from 
DNS for a single mode and verify the quasilinear forms of 
EMFs. 

The above quasilinear EMF terms allow us to identify 
the source of large scale magnetic field growth in a rotating 
plasma with or without radially sheared non-axisymmetric 
perturbations. 


3.2 EMFs and Large Scale Field Growth 

The dominance of the fluctuation-induced quasilinear £z 
in the generation of the large-scale toroidal magnetic field 
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can be seen by examining the mean (averaged in (j> and z) 
toroidal component of the induction equation, ignoring re¬ 
sistivity. This equation is 

^ = -^ + (B-V)VU-(V.V)BU. (21) 

Since there is neither a mean radial magnetic field Br, nor 
velocity field Vr, so the second and third terms on the right 
of Eq. (21) vanish. Note that the second term on the right of 
Eq. (21) is the traditional ”fl effect” which thus vanishes for 
our setup and averaging procedure. The shear (differential 
rotation) does enter through £z, and the first term on the 
right of Eq. (21) is the dominant term. 

Keeping only the first term on the right of Eq. (21) 

= —^§;r) we then see that the three terms on the 
right of Eq. (12) provide distinct paths for large scale fields 
to grow: (1) Radially sheared non-axisymmetric perturba¬ 
tions i. e., the first term in Eq. (12), £zvr', proportional to 
the non-uniformity of the radial displacement of the nonax- 
isymmetric perturbation (2) Uniform non-axisymmetric per¬ 
turbations (stable or unstable) but with background shear 
in the angular flow i. e. the last two terms in Eq. (12). 
We emphasize that all three terms on the right of £z (Eq. 
12) vanish explicitly for axisymmetric modes (m=0 modes, 
Q — 0), as axisymmetric modes are purely growing or de¬ 
caying (i.e. ujr = 0) (Chandrasekar 1961). 

The induction equation for the vertical mean field is 

^ = ^ + (B.V)VU-(V.V)BU. (22) 

Here again, as in equation (21), the last two terms on the 
right vanish and the mean field evolves only through the 
EMF term d£^/dr. For this equation, axisymmetric modes 
can contribute through the radial variations of Vr in £^vr' 
(Eq. 16) and the last term in Eq. (19) to give 

= (y) Szim^O) = 0 (23) 

but only non-axisymmetric modes allow evolution of both 
Bz and B^ through the EMF terms. 


3.3 EMFs appropriate for MRI driven 
fluctuations when kr = Q 

To simplify pinpointing the dominant contributions to the 
EMF from the quasi-linear theory appropriate for MRI 
driven fluctuations for an initially vertical field, we assume 
Vr = 0 (in the limit of kr = 0). This is justified since 
Fig. 3(a) from DNS in the following section shows that the 
term arising from this radial gradient is subdominant. We 
can now solve the quasilinear equations for £r without know¬ 
ing the exact form of the global eigenfunctions. For large 
growth rates (i.e. > uja), Eq. (12) then reduces to 


£zir) 

x( 


'ykBo 


(rn? -\- k^r^) 


^ m?£l'{r)Q.{r)r ^ 2CV,prQ 

TTl o I o 

72 + n 72 + n 


'f- -I- vL{r)‘ 


( 24 ) 
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Figure 3. The profiles of (a) (in black) direct numerical calcu¬ 
lations of total vertical EMF term, LHS of Eq. (12); (in red) 
RHS side of Eq. (12) based on the quasilinear calculations; 
(in blue) the first term in Eq. (12). Linear growth rate of the 
mode at the early phase 'yTori,it ~ 0.9 and (jJrDorbit ~ 2.2 have 

1 f)F 

been used for the calculations of RHS, (b) (B^) and — — S'! 

tiTorbit = 10 and tiTorbit = 10 during the growth of m=l MRI 
mode from DNS with non-zero net flux (solid lines) and zero net 
flux (dashed lines), respectively. The dimensionless magnetic and 
velocity fluctuations at this time of DNS are ~ 0.06 

and {Br/Bol = ~ 0.015, and using DNS values of 

7TA 12, |V^r/VA| 0.015, (Vo/Va = 16). 


This vertical EMF for a single nonaxisymmetic mode and 
the mean flow (e.g. Keplerian) with 7 ~ ilo (where Ho is 
the angular frequency at the inner radial boundary) in the 
kzT > m and 7 ^ > $2 limits are then related by 8z ~ 
-f^‘^\Vrf+2mkBo^L2^\Vr\^ = Q{r)Vr\L Using 

_ ~ 

this equation and B^ ~ “7 7 ^ limit of constant Vr 

and constant mean vertical magnetic field Bq, the large-scale 
toroidal magnetic field can be written as 


B4r) 


rrd Bo 

kzC 

(25) 


1 dQ,{rY 
dlnr 


\V,\^^'SD^[^i^r)fl{r)\'\Vr\ 


showing a direct relationship between differential rotation 
and the generation of large-scale magnetic field. 


4 COMPARING THEORY WITH DNS OF A 
NONAXISYMMETRIC MODE 

From DNS of a nonaxisymmetric mode m = 1, ^^0=12 in 
cylindrical model (Sec. 2), we evaluate (V x 3)2 and com¬ 
pare it to quasilinear calculations of the right side (RHS) of 
Eq. (12) in terms of radial velocity fluctuations Vr = 7^r. 
For the RHS of Eq. (12), the radial velocity fluctuations 
and the eigenvalues from DNS are inserted into the anal^- 
ical forms. Fig. 3(a) shows good agreement between these 
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Figure 4. Radial profiles of time-averaged (b) saturated large- 
scale fields (B^), (B^), (b) J ■ B ( = S (V x B) -(B) from DNS 
of Eq. (1-6) during nonlinear evolution of m=l MRI. 


two calculations. Fig. 3(a) also shows that the first term on 
the RHS of Eq. (12) is subdominant to the last two terms 
which depend on the mean flow. 

We have verified the dominance of the first term on the 
right of Eq. (21) from DNS of a single-mode m=l MRI. The 
large-scale starts to grow, even when initially zero, as the 
instability develops. Figure 3(b) shows {B^) as computed 
from the DNS during the linear phase of single-mode simula¬ 
tions with non-zero net flux and the first term on the RHS of 
Eq. (21), right before the saturation, as also measured from 
the DNS. As seen, the mean toroidal field is correlated with, 
and directly generated by the vertical EMF. Similarly, the 
mean generated in the net-zero flux simulations shown 
in Fig. 1 is also correlated with the vertical EMF. Fig 3(a) 
shows that the main contribution to the EMF comes from 
the last two terms of Eq. (12). Thus B^ ~ ^ 

is directly dependent on the shear-flow (V^) or differential 
rotation {d£l/dr) in the presence of a finite amplitude fluc¬ 
tuation. 

In addition to the EMF components themselves, the 
magnetic field-aligned EMF plays important role in the 
quasi-linear regime as the mode starts to saturate. Evo¬ 
lution of both components (toroidal and vertical) of the 
large-scale magnetic field is only possible for nonaxisymmet- 
ric fluctuations because only in this case are both and 
£ri> (Eq. 12, 16) non-zero. The generation of the large-scale 
toroidal field is due to the vertical EMF and the redistri¬ 
bution (and amplification around r=0.2-0.4) of Bz is due 
to the nonzero toroidal EMF. Figure 4 shows the profiles 
of time-averaged saturated large-scale toroidal and vertical 
fields as well as the EMF parallel to B. 

Simulations with axisymmetric fluctuations also show 
the amplification of vertical large-scale field (Bz) but with¬ 
out a generation of B^. The amplification of Bz from ax¬ 
isymmetric modes (Eq. 23), which may contribute to the 



Figure 5. Field line visualizations during nonlinear MRI simula¬ 
tions in a cylinder (a) 2-D (b) 3-D when m=l mode perturbation 
is dominant (c) toroidal top view in 3-D shows the twisting of the 
field lines. 

mode saturation in a cylinder (Ebrahimi et al. 2009) results 
from the curvature terms and is absent for channel modes 
in a local Cartesian model, as discussed below. 


5 VISUALIZATIONS FROM NUMERICAL 
SIMULATIONS 

The physical picture of generating (B^) can further be exam¬ 
ined through comparing visualizations of the field lines with¬ 
out and with nonaxisymmetric MRI perturbations in 2-D 
and 3-D as shown in Fig. 5. For toroidal m=0 perturbations, 
weak vertical magnetic field lines are toroidally stretched 
(Fig. 5(a)) according to the second term in Eq. (21). How¬ 
ever due to toroidal symmetry, (B^) = 0 as the positive 
and negative contributions from the perturbations remain 
on the same vertical surface, and only the mean vertical 
field is amplified through dS^jdr. (Ebrahimi et al. 2009). 
In the presence of nonaxisymmetric perturbations in 3-D 
nonlinear simulations however, the field lines are stretched 
and twisted (Fig. 5 b,c). As a consequence, {B^j,) yf 0 since 
now the oppositely signed toroidal field contributions from 
perturbations are displaced radially from one another. Since 
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{Br) is zero, the standard “f2 effect” [{Br)V^ contribution 
in Eq. (21) is zero. 
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6 LOCAL CARTESIAN QUASILINEAR 
ANALYTIC CALCULATIONS OE EMFS 


6.1 General form of the EMEs 


Here, we present the analog to Eqs. (7-12) in local Carte¬ 
sian coordinates (x, y, z) in a frame rotating with fixed an- 
gnlar velocity fl = floez and a linear shear velocity of 
Vo = Vy{x)ey. We again assume a vertical field Bo~z but now 
assume perturbed velocity and magnetic field of the form 
i{x,y,z,t) = \^^{x),e,y{x),i^{x)]exp(;^^t-ikyy + ik^z), 

where, 7c = 'y + iuJr- In this rotating unstratified system, we 
include the Coriolis force and the centrifugal force, and again 
assume that in equilibrium, the latter is canceled by the ra¬ 
dial pressure gradient. (As for the cylindrical case, the role 
of gravity vs. radial pressure gradient are interchangeable 
for our incompressible, unstratified case. ) The momentum 
equation = —pV.VV-|-2(nx V)+fIxnxr-|-J xB—VP 
and the induction equation 9tB = V x (V x B) are linearized 
in the incompressible limit to give: 


7^7c 2HoV)/ — ikzBoBx / p ^ 

yVy + Vy{x)Vx + 2Q.oVx = ik^BoBy/ p + ikyX 

7I4 = ikzBoBz/p — ikzX 


(26) 


and 

'yBx — 'ikzBoVx 

^By = Vy{x)'Bx -\- ikzBoVy (27) 


'^Bz = ikzBoVz, 


where primes denote variation in x direction, {e.g. d/dx) 
Using Vx = 7 ^ 1 , Vy = yCy-VyOC and 14 = Eqs. (26) 
and (27) can be written in terms of the displacement vector 
as 

o o ' dX 

[ 7 ^ +iDA + 2rQ,oVy (a:)]^!, - 2'yQ,o^y = - — (28) 

( 7 ^ -I- + 2^Q.o^x = ikyX (29) 

( 7 ^ -I- oj\)^z = -ikzX. (30) 

Analogous to the cylindrical case (Sec. 3), from these 
sets of equations, the quasilinear vertical EMF in terms of 
the Eulerian velocity fluctuations (V = d^/dt -|- V ■ ^ ■ 

VV), 14 = jC is reduced to: 


£z — (V X B) 2 , — £ ZVr' + £Z(shear) + f ZH, (31) 

where 




ky'fkzBo 

{kl + kl) 


1414* 


7^ -I- i9^(x) 


(32) 


Z{shear) 


ky^kzBo 

■ (fc2 + kl) 


kyVy{x) {kyVyix) U!r) 


(72 _|_ 


114 




and 


£zno = 


2'ykxBo£lo'&{x), 2 


G2 


(7 -|-d(a:) -uja)\Vx 




(33) 

(34) 


where 'd{x) = kyVy{x) — ojr, 7 = 7 — ii9(x), and = (ky + 
kz){C + C{x))[4:'y'^g'^{x) + ( 7 ^ -I- y^(x) + CaT]- 2ls seen in 
these equations, in the local Cartesian model, the rotation 
and shear are independent. 

The first contribution £ zvx' shows how nonaxisym- 
metric perturbations with radial shear, even without any 
explicit mean shear flow can source a vertical EMF. In 
the absence of rotation, the second contribution, £z(shear), 
shows a direct dependence of vertical EMF on the lin¬ 
ear shear. A mean shear flow, Vy{xY combined with a 
radially uniform non-axisymmertic {ky 7 ^ 0 ) perturbation 
is sufficient to produce £z(shear)- The last contribution 
£zn, which vanishes in the absence of angular velocity, 
shows that a contribution to the vertical EMF can result 
from a finite angular velocity (Ho) for non-axisymmetric 
perturbations. 

Similarly the azimuthal EMF is given by 


£y ^ 13 )y £ yVr' T £y{shear) T £yCLQi 

where 

(1414*') 

(72 _|_ i[}'i{x)) ’ 


£ yv x’ 'yBo 


1 - 


kt 


{kl + kl) 


(35) 


(36) 


£ y(shear) 


- 7 B 0 



Cy 

{kl + kl) 


kyVy{xy{kyVy{x) - UJr)\Vx\^ 

(72 _|_ 

(37) 


and 

= 27fc,^-SoHoMfc.l4(2:)-02.) ^2 _ ^2 _ ^2(^))|^|2_ 

(38) 

Spatial derivatives of Eq. (35) could, in principle, also 
generate and redistribute the vertical field Bz due to non- 
axisymmetric modes {ky 7 ^ 0). However, the fastest growing 
axisymmetric modes-the channel modes {ky = fcj, = 0 )— 
do NOT contribute in either the vertical or azimuthal EMF 
(Eq. 31 and Eq. 35) obtained above. In contrast, for the 
global cylindrical model, even for radially uniform axisym¬ 
metric modes {kr = m = 0), the last term in Eq. (23) DOES 
contribute in the amplification of vertical field. This dis¬ 
tinction highlights at least one circumstance in which the 
absence of curvature in the Cartesian model removes a con¬ 
tribution that could be present in the global rotator. 


6.2 Exact expression for large scale field with 
only linear shear when kx = 0 


A large-scale magnetic field can be generated by the EMF of 
Eq. (31) from any of the independent contributions in (32- 
34). In the absence of rotation, a large-scale magnetic field. 

By ~ ~ 7 , can directly be generated via a linear flow- 

shear and a radially uniform non-axisymmetric {ky,kz Y Q, 
kx = 0 ) perturbation. 


By{x) = 


{kl + kl) 


kyVy{x)' {kyVy{x) — Ur 


(72 _j_ 


IK 


(39) 
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Figure 6. generated with positive shear flow {V^{r)/VA = 
80(r/a)®) and large initial amplitude forcing (solid line); and only 
with forcing amplitude varying with radius (dashed line). 

This is an exact analytical equation for a large-scale az¬ 
imuthal magnetic held generated via a linear mean shear- 
flow and any perturbations with nonzero ky and kz- 

The large-scale field given in Eq. (39) is consistent with 
previous studies of large scale field growth from the combina¬ 
tion of linear shear with randomly forced turbulence (Vish¬ 
niac & Brandenburg 1997; Yousef & et al. 2008; Heinemann 
et al. 2011; Mitra & Brandenburg 2012; Sridhar & Singh 
2014). However, our calculations explicitly reveal the most 
minimalist conditions needed for growth in the absence of 
rotation: a background linear shear and an imposed non- 
axisymmetric perturbation with nonzero ky,kz. Helical ve¬ 
locity perturbations are not required. 

Generation of By in this case of mean shear can be vi¬ 
sualized by considering an perturbation in the x direction 
and then considering why both z and y variations are needed 
to produce a net field in vertically averaged planes. If there 
were no z variation in the perturbation then the mean shear 
would produce no toroidal field even before vertically aver¬ 
aging. And if there were a vertical variation but no y vari¬ 
ation, then the the mean shear would produce By{x) = 0 
from vertical averaging. 


7 LARGE-SCALE FIELD IN THE CASE OF 

STABLE FLOW 

Our quasi-linear theory imposes fluctuations and back¬ 
ground shear as a starting point whereas in DNS, the fluc¬ 
tuations can directly result from the MRI. The quasilinear 
theory shows that the growth of B^ in Eq. (21) via drSz 
does not requires a shear profile favorable to the MRI, just a 
source of fluctuations and differential rotation of either sign. 
To show this, we numerically computed the large-scale field 
growth from quasilinear theory by initializing single mode 
fluctuations (/m,fe(r, 0) with a polynomial dependence on r) 
in the simulations and forcing amplitudes of 10% on top of 
a stable equilibrium flow dfl/dr > 0. 

As the fluctuation energy decays, B^ grows via Ez from 
(Eq. 12). Figure 6 shows the large-scale toroidal field gen¬ 
erated using V 4 ,{r)IVA = 80(r/a)®. The profiles are time- 
averaged during the decay phase. In Figure 6 we have also 
shown the case when B^ is generated by forcing only with 
the same radially dependent fluctuations but in the absence 
of mean shear. For this latter case, £z in Eq. (12) is then 
dominated by the first term on the right. Comparing the 


two cases, we see that for small r/a the case with only 
radius dependent fluctuations (dashed line in Fig. 6) also 
captures the growth of B^ as for the case with both fluctua¬ 
tions and shear. But for radii of large shear rja > 0.75, the 
flow-dependent EMF terms (the last two terms of Eq. 12) 
dominate. 


8 SUMMARY AND CONCLUSIONS 

In summary, we have shown from both numerical simula¬ 
tions and semi-analytic quasi-linear theory how radially al¬ 
ternating large-scale toroidal fields averaged vertically and 
azimuthally can be generated from MHD flow-driven fluctu¬ 
ations. These fields are found in MHD DNS for both zero- 
net-flux and non-zero-net-flux initial configurations in both 
the quasi-linear regime and the fully saturated non-linear 
regime. 

Given nonaxisymmetric fluctuations (with nonzero ver¬ 
tical and azimuthal perturbations), we calculated the con¬ 
tributions to the quasilinear fluctuation-induced EMFs in 
both cylindrical and Cartesian coordinates. We have sepa¬ 
rated the derivation of the global and local models so that 
the reader can study them separately. 

We have not presented physical interpretations of all cir¬ 
cumstances that can lead to growth from these equations, 
but have provided the general forms of the EMFs and iden¬ 
tified the minimum requirements for growth. Table 1 sum¬ 
marizes these requirements for a nonzero EMF and large- 
scale field growth in both cylindrical and Cartesian models. 
In general, we find a direct relationship between dynamo 
generating EMFs and differential rotation in the cylindri¬ 
cal model, or linear shear in the local Cartesian model. The 
vertical EMF associated with fluctuations in the presence of 
an initial vertical field is sufficient to generate an azimuthal 
large-scale field for non-zero differential rotation (in rotat¬ 
ing system) but requires non-zero flow shear in the local 
Cartesian model for a non-rotating system. 

Table 1 also highlights that due to the absence of cur¬ 
vature terms, the local Cartesian model is more restrictive 
for field growth than global cylindrical model. According to 
our Cartesian EMF calculations, the fastest growing channel 
modes (with ky — k^ = 0) (Goodman & Xu 1994) found in 
shearing box simulations do NOT contribute to the EMFs 
in the local approximation (and thus the saturation of these 
modes) but the analogous modes can amplify large-scale 
fields and contribute to MRI saturation in global cylindrical 
simulations (Ebrahimi et al. 2009). 

In the case of a large scale flow-driven instability, the 
free energy source from the large scale motion can be the 
source of the needed fluctuations. For the global cylinder, 
we have indeed found explicit dynamo generation of {B^) 
from DNS where the MRI produces a fluctuation-induced 
vertical EMF Ez- The DNS provide properties of the fluctu¬ 
ations that we use as inputs to a quasi-linear calculation of 
the dynamo growth for a single mode. The DNS large scale 
field growth and the associated quasi-linear dynamo calcu¬ 
lations are in reasonable agreement. Our study of the single 
mode evolution its correspondence with DNS highlights that 
’’turbulence” (defined as non-linear mode coupling) is not 
actually essential for the large-scale field growth and that 
insight is gained even from single mode analyses. 
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EMF radial derivative 

Restriction Maintains Finiteness? 

global cylinder (with > 0 ; 7 ^ 0 ; /c 7 ^ 0 ; fcr 7 ^ 0 ; m > 0 ) 

kr = 0 

n'{r) = 0 

m = 0 

e'^{T) + 0 

yes 

yes 

no 


yes 

yes 

yes 

£'^{r) and ^’^(r') 7 ^ 0 (dynamo) 

yes 

yes 

no 

local Cartesian (with > 0; 7 ^ 0; /c 7 ^ 0; /cx 7 ^ 0; ky > 0) 

kx =0 

II 

0 

ky = 0 


yes 

no 

no 

e'y{x) / 0 

yes 

no 

no 

£y{x) and 7 ^ 0 (dynamo) 

yes 

no 

no 

local Cartesian (with Qq = 0; 7 ^ 0; /c 7 ^ 0; 7^ 0; ky > 0) 

kx = 0 

V^ix) = 0 

ky = 0 


yes 

no 

no 

£'y{x) 7 ^ 0 

yes 

no 

no 

£y{x) and £'^(x) 7 ^ 0 (dynamo) 

yes 

no 

no 


Table 1. Table summarizing the minimal ingredients needed to maintain the finiteness of the two components of the EMFs separately 
and together. The first column indicates three general cases: the global cylinder and two Cartesian cases. To read the table for each of 
these cases, consider the global cylinder case as an example: the second row of the first column indicates general ingredients that our 
global cylinder could have (rotation > 0, differential rotation fl' yf 0, and general perturbations with A: yf 0, fcr yf 0, and m > 0 for 
non-axisymmetry). Each of the subsequent columns indicate a restriction that reduces this generality. The entries ”yes” or ”no” within 
these columns indicate whether the corresponding quantity in the first column is finite when that restriction is imposed. For each of the 
three cases (global cylinder, Cartesian with flo > 0 and Cartesian with Qq = 0) listed, a first and second row provide the information on 
the finiteness of each component of the and the third row provides the information on the finiteness of both components together. This 
latter circumstances is needed to supply the large scale dynamo. Note that if either of the two separate components has a ”no” entry 
then that also implies that both together cannot be finite. The conditions for finiteness (and thus dynamo action) are different in the 
two geometries because curvature terms are absent in the Cartesian approximation. 


Our results also show that the traditional “fl effect” of 
shear on the mean field is absent when the initial mean field 
is vertical and the averaging is over vertical surfaces. Instead 
the essential shear operates on the fluctuations. The field 
growth can be entirely described by working with the EMF 
directly, non-axisymmetric (though not necessarily helical) 
velocity perturbations are essential for large scale growth as 
evidenced from direct visualization of the field lines in DNS 
and from the quasi-linear theory. We should also note that 
in much of the MRI dynamo literature, large-scale fields in 
shearing boxes are computed via planar averages (and av¬ 
eraged over the direction of the nonuniformaty of the mean 
flow) leaving mean fields as a function of z direction. There, 
because of the averaging and boundary conditions for the 
shear box “fl effect” can still survive. Here, our averaging is 
over vertical surfaces, and not along the direction of mean 
flow variation. An important lesson is that the averaging 
procedure and boundary conditions have important impli¬ 
cations for the dominant contributions to the EMF. 

By calculating the complete form of EMF for both 
global cylindrical and Cartesian cases, we have demon¬ 
strated the minimum ingredients for large scale field growth 
in both of these two models. Our results suggest that the 
quasilinear and nonlinear fluctuation-induced EMF may 
provide fundamental insight into the growth and sustenance 
of large-scale dynamo in these flow-driven systems. The cal¬ 
culations herein provide a more general approach to iden¬ 
tifying the origin and minimal ingredients needed for large 
scale dynamo growth in unstratified rotating and differen¬ 
tially rotating systems or linearly sheared systems. 

Although we leave a detailed analysis making explicit 
connections to previous approaches of incoherent alpha 
(Vishniac & Brandenburg 1997; Brandenburg 2005; Bran¬ 
denburg et al. 2008; Mitra & Brandenburg 2012; Sridhar & 
Singh 2014) and or shear current effects as an opportunity 


for further work, we emphasize two points in this context. 
First we have intentionally avoided using the a formalism 
and worked directly with only the EMF. Second, we find 
that the absolute minimum conditions for radially depen¬ 
dent large scale field growth are non-axisymmetric velocity 
fluctuations plus linear shear. The velocity fluctuations do 
not need to be helical at any time. In this way our global 
and local calculations provide a more minimalist set of con¬ 
ditions for growth than the that of a fluctuating kinetic helic- 
ity (Vishniac & Brandenburg 1997). We note however, that 
in the quasi-linear regime, the large scale magnetic field does 
(as a function of radius) (see Fig. 4) develop a field aligned 
EMF, which is a source term for sum of the time derivative 
of large scale magnetic helicity and divergence of large scale 
helicity flux, as previously confirmed for the global cylin¬ 
drical case (Ebrahimi & Bhattacharjee 2014). Here we have 
not studied the non-linear/saturating effects of the growth 
of small scale magnetic helicity helical fluctuations, nor the 
EMF and mean magnetic field correspondence during the 
nonlinear saturation. More detailed calculations for the non¬ 
linear phase of DNS (by P. Bhat et.al in preparation) do 
show a direct correlation of large-scale field with the EMF 
terms in the nonlinear regime that we have presently com¬ 
puted only in the quasli-linear approximation. 

Finally, we note that our large scale fields show radial 
reversals and these would be sites of current sheets. If we 
think toward generalizations to stratified rotators that form 
coronae, only magnetic structures of large enough scale sur¬ 
vive buoyant rise into coronae where they can dissipate and 
transport angular momentum non-locally (Blackman & Pes- 
sah 2009). If our present toroidal field structures and rever¬ 
sal scales survive stratified generalizations, they provide a 
scale for coronal structures and current sheets that link the 
large scale field directly to structures associated with coro¬ 
nal transport and dissipation. 
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